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ABSTRACT 

We present results from a wlOO ks Chandra observation of the 2QZ Cluster 1004+00 structure at z = 2.23 
(hereafter, 2QZ Clus). 2QZ Clus was originally identified as an overdensity of four optically-selected QSOs at 
z = 2.23 within a 15 x 15 arcmin^ region. Narrow-band imaging in the near-IR (within the K band) revealed 
that the structure contains an additional overdensity of 22 z = 2.23 Ha-emitting galaxies (HAEs), resulting 
in 23 unique z = 2.23 HAEs/QSOs (22 within the Chandra field of view). Our Chandra observations re- 
veal that 3 HAEs in addition to the 4 QSOs harbor powerfully accreting supermassive black holes (SMBHs), 
with 2-10 keV luminosities of w(8-60) x lO"^ ergs s"' and X-ray spectral slopes consistent with unobscured 
AGN. Using a large comparison sample of 210 z = 2.23 HAEs in the C/ianc/ra-COSMOS field (C-COSMOS), 
we find suggestive evidence that the AGN fraction increases with local HAE galaxy density. The 2QZ Clus 
HAEs reside in a moderately overdense environment (a factor of ss2 times over the field), and after exclud- 
ing optically-selected QSOs, we find the AGN fraction is a factor of sa3.5^2 2 times higher than C-COSMOS 
HAEs in similar environments. Using stacking analyses of the Chandra data and Herschel SPIRE observations 
at 250/im, we respectively estimate mean SMBH accretion rates (Mbh) and star-formation rates (SFRs) for 
the 2QZ Clus and C-COSMOS samples. We find that the mean 2QZ Clus HAE stacked X-ray luminosity is 
QSO-like (L2-10 keV « [6-10] x lO"^ ergs s"'), and the implied Mbh/SFR w (1.6-3.2) x 10"^ is broadly consis- 
tent with the local Mbh/M,^ relation and z~ 2 X-ray selected AGN. In contrast, the C-COSMOS HAEs are on 
average an order of magnitude less X-ray luminous and have Mbh/SFR w (0.2-0.4) x 10"^^, somewhat lower 
than the local Mbh/M* relation, but comparable to that found for z ~ 1-2 star-forming galaxies with similar 
mean X-ray luminosities. We estimate that a periodic QSO phase with duty cycle w2-8% would be sufficient 
to bring star-forming galaxies onto the local M^u/M^ relation. This duty cycle is broadly consistent with the 
observed C-COSMOS HAE AGN fraction (wO.4-2.3%) for powerful AGN with Lx ^ lO"'* ergs s"'. Future 
observations of 2QZ Clus will be needed to identify key factors responsible for driving the mutual growth of 
the SMBHs and galaxies. 

Subject headings: cosmology: observations — early universe — galaxies: active — galaxies: clusters: general 
— surveys — X-rays: general 



1. INTRODUCTION 

Successful theoretical models characterizing the formation 
and growth history of galaxies and supermassive black holes 
(SMBHs) generally require that feedback from active galactic 
nuclei (AGN) play a crucial role in regulating the growth of 
galaxy bulges (e.g., De Lucia etal. 2005; Bower etal. 2006, 
2008; Croton et al. 2006; Fanidakis et al. 2012), leading to the 
local M-cr and Mbh-Mj, relations (e.g., Bennert etal. 2011). 
AGN feedback processes are predicted to be most important 
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in the most massive elliptical galaxies formed within the most 
massive dark-matter halos. However, in the local universe, 
these same elliptical galaxies contain the oldest stellar popu- 
lations and the most dormant SMBHs in terms of their spe- 
cific mass accretion rates. Mass accretion in these systems 
is likely limited by the presence of a hot interstellar medium 
and periodic mechanical feedback from radio powerful AGN 
(e.g.. Best etal. 2005; McNamara & Nulsen 2012). There- 
fore, directly observing the growth of these systems requires 
challenging observations of very distant precursors to today's 
ellipticals. From models of large-scale structure formation, it 
is predicted that growth of modem ellipticals takes place in 
high-density regions at z 2-3 (e.g., Kauffmann etal. 1996; 
Governato etal. 1998; Volonteri etal. 2003; De Lucia etal. 
2006). 

Thus far, efforts to observe the distant z ~ 2-3 galaxy proto- 
clusters SSA22 and HS 1700+64, likely pre-cursors to rich lo- 
cal clusters like Coma (e.g., Steidel et al. 1998, 2005), have re- 
vealed that the AGN fraction is a factor of w2-16 times larger 
(la statistical range) in the protocluster environment com- 
pared with the low-density field (e.g., Lehmer etal. 2009a,b; 
Digby-North etal. 2010). This contrasts with local galaxy 
clusters, which contain lower AGN fractions in their cores 
(e.g.. Martini etal. 2009; Ehlert etal. 2012). The enhanced 
SMBH growth rate in the protocluster environment is likely 
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Fig. 1 . — Evidence for wide-spread AGN activity for z = 2.23 HAEs in the 2QZ Clus structure. We show the adaptively-smoothed 0.5-7 keV Chandra image 
of the 2QZ Clus field. For illustrative purposes, the image was binned by a factor of 4 (in both RA and DEC) and has been exposure-corrected by dividing the 
smoothed image by a smoothed exposure map (see §2.2). The aim point of the image has been indicated with a cross. The z = 2.216-2.248 QSOs that were used 
to select the 2QZ Clus field have been highlighted with green squares (Croom etal. 2001, 2004). The WFCAM HAE survey region (dashed square) and z = 2.23 
source candidates {red circles) have been highlighted (see Matsuda etal. 2011). BzK candidates have been highighted with blue diamonds; these sources have 
colors consistent with actively star-forming galaxies at z 1.4—2.5. 



linked to the presence of larger gas reservoirs being fed onto 
more massive galaxies (e.g., Steidel et al. 2005) and SMBHs, 
as well as an increase in mergers before the cluster becomes 
virialized. These structures are therefore ideal for observing 
and characterizing the interdependent processes involved in 
the growth of the most massive SMBHs and galaxies, and 
follow-up studies of protocluster AGN have revealed inter- 
esting first results. For example, in the SSA22 protoclus- 
ter at z = 3.09, many of the AGN are found to be coinci- 
dent with large-scale (^100 kpc) Lyman-a emitting nebulae 
(Geach et al. 2009) that may be indicative of emission-line gas 
being powered by AGN. 

To expand upon our efforts to understand AGN activity in 
the high-density environment at high redshift, we have per- 
formed Chandra observations of an overdensity of actively 
star-forming galaxies and optically-selected QSOs at z = 2.23: 
the 2QZ Cluster 1004H-00 structure (hereafter, 2QZ Clus). 
2QZ Clus was discovered by Matsuda etal. (2011) by first 
searching the 2dF QSO Redshift Survey (2QZ; Croom et al. 
2001, 2004) and identifying an overdensity of quasars within 
the redshift interval z = 2.216-2.248, corresponding to the Ha 
line in the UKIRT WFCAM H2SI filter, and then performing 
narrow and broad band near-infrared imaging (via the H2SI 
Af = 2.121 ^va and /T-band filters) to identify star-forming ac- 
tive galaxies at z = 2.23 via the Ha emission line. The redshift 
range of the Ha line corresponding to the 50% transmission 
wavelengths of the H2SI band is z = 2.216-2.248. These Ha 



emitters (HAEs) therefore have both redshift identifications 
and Ha-based measures of their star-formation rates (SFRs). 
Matsuda et al. (2011) revealed that the 2QZ Clus field con- 
tains an overdensity of 22 vigorous star-forming HAEs (Ha- 
based SFR ^ 14 Mq yr~' without any correction for ex- 
tinction) in a 13.7 X 13.7 arcmin^ region («22 x 22 comov- 
ing Mpc^). These observations probe the tip of the iceberg 
of the actively star-forming z = 2.23 2QZ Clus galaxy popula- 
tion with additional contributions from luminous AGN. Iden- 
tification of the lower-SFR and Ha-obscured population is 
still needed to fully characterize the large-scale nature of this 
structure. 

In this paper, we present first results from the Chandra 
observations of 2QZ Clus. In § 2, we discuss our analysis 
methods and provide reduced data products, including im- 
ages, exposure maps, and point-source catalogs. In § 3, we 
present the properties of the X-ray detected sources in the 
2QZ Clus structure, and compare these with an equivalent 
unbiased sample of z = 2.23 HAEs found in the COSMOS 
survey field from the High-z Emission Line Survey (HiZELS; 
Geach et al. 2008; Sobral et al. 2009). In § 4, we compare the 
relative growth rates of 2QZ Clus SMBHs and galaxies. We 
further discuss the nature of the 2QZ Clus structure itself and 
put the AGN activity into the context of the broader z = 2.23 
population of HAEs found in COSMOS. Finally, in § 5, we 
summarize our results. 
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TABLE 1 

Chandra POINT-SOURCE CATALOG OF THE 2QZ Clus Field 



Position (J2000) 



Net Counts 



Exposure (ks) 



Count-Rate (10~* cuts s"') 



Source ID aj2000 <5j2ooo 

(1) (2) (3) 

1 10 03 00.7 +00 10 54.1 

2 10 03 00.8 +00 11 25.1 

3 10 03 05.9 +00 09 49.2 

4 1003 09.2 +00 11 31.1 

5 10 03 09.3 +00 12 24.3 

6 10 03 09.6 +00 09 01.1 

7 10 03 09.8 +00 13 12.9 

8 10 03 13.1 +00 13 05.1 

9 10 03 13.2 +00 10 26.8 

10 10 03 13.9 +00 12 17.5 



FB 
(4H6) 



SB 

(7)-(9) 



HE 

(10H12) 



FB 

(13) 



SB 
(14) 



HB 
(15) 



FB 
(16)-(18) 



SB 

(19)-(21) 



HB 

(22)-(24) 



I (: 4+10.0 

-8.9 
4.7 7+14.2 
^^■^-11.4 
67 3+'2 ' 

24.7:8-9 

64 2+"-3 
°^-^-10.2 

68.9:1':^ 
14.01:° 

C7 7+10.2 
-"■'-9.1 

41-91:? 
13.7;^:° 



<20.2 

^'-■^-1.1 

<15.5 

i4.4;f| 

27.8!^:5 
54.5_g:3 
U3^l 
29.0:^:? 
21.6+«9 



11.3^ 



i2_ 



<26.0 
<36.3 

66.8:li:^ 

<22.4 
38.3!94 
<21.8 
<19.7 
30.4:|l 

2o.l1:? 

<15.6 



53.5 
24.6 
61.7 
71.6 
69.5 
46.0 
65.9 
74.5 
76.3 
78.7 



47.7 
21.9 
57.1 
68.7 
66.1 
42.5 
61.8 
71.9 
74.1 
77.0 



77.2 
35.8 
80.5 
83.3 
83.0 
60.3 
82.6 
84.6 
85.1 
85.9 



9 9+2.0 
^■^-1.8 

11 3+2' 
3 4+'-3 
9.2!l'^ 



-1.5 
I5.1II 

n+1.2 
-1.0 
7 q+1.4 
'■^-1.3 
5.5+1-3 



2.0+ 



1.7 



-1.1 
+0.9 
-0.8 



<4.5 
13.2+5-2 

<2.8 
T 1+0.9 
^-'-0.8 
4 2+' -2 

1 3 2+2-2 

T q+1.0 
^-^-0.8 
4 l+'O 

2 Q+0-9 



1.5 



<3.7 
<10.8 
8.8!|-^ 

<2.8 
4 6+' 2 

<3.6 

<2.4 

3 6+'-° 
^-"-0.8 

7 4+0-9 
<1.9 



NOTE. — [A portion of Table 1 is shown here to illustrate content. All 133 rows and 28 columns are provided in the electronic version.} Col.(l): Chandra source 
ID for 2QZ Clus field. Col. (2) and (3): Right ascension («j20oo) and declination (5j2ooo), respectively. Right ascension is quoted in units of degrees, minutes, 
seconds. Declination is quoted in degrees, arcminutes, and arcseconds. Col.(4)-(12): Net counts {N) and Ict upper/lower bounds, computed following §2.2, 
for FB, SB, and HB. Col.(13)-(15): Vignetting-corrected effective exposure times (in kiloseconds) for FB, SB, and HB, respectively. Col.(16)-(24): Vignetting 
corrected count-rates and \ a upper/lower bounds, computed following the methods described in §2.2, for FB, SB, and HB. Col.(25)-(27): X-ray flux in units of 
ergs cm"2 g-l for the FB, SB, and HB, respectively. Col. (28): Notes related to source ("M" = Manual photometry required due to image edge; "H" = Coincident 
with z = 2.23 HAE; "Q" = Coincident with optically identified QSO). 



The Galactic column densities for 2QZ Clus and COSMOS 
are 3.0 x 10^" cm"^ and 2.5 x 10^" cm"^, respectively. All of 
the X-ray fluxes and luminosities quoted throughout this pa- 
per have been corrected for Galactic absorption. In the X-ray 
band, we make use of three bandpasses: 0.5-2 keV (soft band 
[SB]), 2-7 keV (hard band [HB]), and 0.5-7 keV (full band 
[FB]). Values of //,) = 70 km s"' Mpc"', = 0.3, and 
= 0.7 are adopted throughout this paper (e.g., Spergel etal. 
2003). 

2. OBSERVATIONS, REDUCTION, AND ANALYSIS 

Throughout this paper we will be comparing the SMBH 
growth among HAEs in the 2QZ Clus structure with HAEs 
in the COSMOS field. To do this effectively, we designed our 
Chandra observations of 2QZ Clus to be of comparable depth 
and quality as the Chandra data products that are akeady 
available from the Chandra COSMOS survey (C-COSMOS; 
Elvis etal. 2009; Puccetti etal. 2009; Civano etal. 2012). 
Therefore the Chandra analysis that follows has been cus- 
tomized to produce catalogs and data products that can be 
directly compared with those available from C-COSMOS. 

2.1. Data Reduction 

We obtained a wlOO ks Chandra exposure consisting 
of a single 16.9' x 16.9' ACIS-I pointing (Chandra Obs- 
ID 13976; taken over the 13th and 14th of January 2012; PI: 
B.D. Lehmer) centered on the 2QZ Clus region surveyed by 
Matsuda etal. (201 1; see Fig. The observation was cen- 
tered on the aim point coordinates aj2ooo = 10:03:43.3 and 
<5j200() = +00:13:26.47 and was oriented at a roll angle that 
was 56.9 deg from north. The total duration of the obser- 
vation was 99.6 ks. For our data reductions, we made use of 
CIAO V. 4.4 with CALDB v4.5.0. We began by reprocessing 
our events lists, bringing level 1 to level 2 using the script 
chandra_repro. The chandra_repro script runs a 
variety of CIAO tools that identify and remove events from 

'" Note that during the obsei-vation, the ACIS-S3 chip was on; however, 
due to its large off-axis angle and non-coincidence with any z ~ 2.23 sources 
in the 2QZ Clus field, we chose to exclude data from the ACIS-S3 CCD. 



bad pixels and columns, and filter the events list to include 
only good time intervals without significant flares and non- 
cosmic ray events corresponding to the standard ASCA grade 
set (ASCA grades 0, 2, 3, 4, 6). 

Using the reprocessed level 2 events list, we constructed a 
first FB image and a point-spread function (PSF) map (using 
the tool mkpsfmap), which corresponded to a monochro- 
matic energy at 1.497 keV and an encircled counts fraction 
(ECF) set to 0.393. We constructed an initial source cata- 
log by searching our FB image with wavdetect (run with 
our PSF map), which was set at a conservative false-positive 
probability threshold of 1 x 10"^ and run over seven scales 
from 1-8 (using a \/2 sequence: 1, V2, 2, 2V2, 4, 4\/2, 
and 8). To sensitively measure whether any significant flares 
remained in our observations, we constructed point-source - 
excluded 0.5-7 keV background light curves for the observa- 
tion in a variety of time bins (spanning 10-800 s bins). We 
found no evidence of any significant ( ^ 5 ct) flaring events 
throughout the observation, and therefore considered our re- 
processed level 2 events file to be sufficiently cleaned. 

Next, using the initial X-ray source catalog and an optical 
I < 22 mag source catalog from the Sloan Digital Sky Sur- 
vey Data Release 6 (SDSS-DR6; Adelman-McCarthy etal. 
2008), we registered our aspect solution and events list to the 
SDSS-DR6 frame using CIAO tools repro ject_aspect 
and repro ject_events, respectively. The resulting as- 
trometric reprojections gave very small astrometric adjust- 
ments, including linear translations of 6x = -Q.l2 pixels and 
Sy = +0.30 pixels, a rotation of -0.0097 deg, and a pixel scale 
stretch factor of 1.00018. 

2.2. Point-Source Catalog Production 

Using the reprojected aspect solution and events file (see 
§2.1), we constructed standard-band images, 90% ECF PSF 
maps (corresponding to 1.497 keV, 4.51 keV, and 2.3 keV 
for the SB, HB, and FB, respectively), and exposure maps 
in the three standard bands. The exposure maps were con- 
structed following the basic procedure outlined in § 3.2 of 
Hornschemeier etal. (2001); these maps were normalized to 
the effective exposures of sources located at the aim points. 
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Fig. 2. — Observed X-ray spectral shapes (HB-to-SB count-rate ratio) 
vs. 2-10 keV luminosity for X-ray detected z = 2.23 HAEs and QSOs in 
2QZ Clus and C-COSMOS fields. The red open circles and green open 
squares show 2QZ Clus HAEs and z = 2.23 QSOs, respectively; blue di- 
amonds indicate sources that satisfy the BzK star-forming galaxy criteiia. 
C-COSMOS z = 2.23 HAEs are shown as filled circles. Results from stacking 
HAE samples in different overdensity bins in C-COSMOS are shown as black 
circles with the letter "C" and stacking of 2QZ Clus HAE samples including 
and excluding pre-selected QSOs used to identify the 2QZ Clus structure are 
respectively shown as black and gray squares with the letter "Q" (see §4.2 
for details). The hoiizontal dashed line shows the expected count-rate ratio 
for an unobscured AGN with a power-law SED of F = 1.9. The shaded re- 
gion shows the expected count-rate ratio for AGN obscured by a column of 
log(Nii/cnr-) ^ 23.5. All of the 2QZ Clus sources have spectral shapes con- 
sistent with luminous unobscured AGN, while the C-COSMOS AGN span a 
broader range of spectral slopes. 

This procedure takes into account the effects of vignetting, 
gaps between the CCDs, bad column and pixel filtering, and 
the spatially dependent degradation of the ACIS optical block- 
ing filter. A photon index of F = 1 .4 was assumed in creat- 
ing the exposure maps. In Figure 1, we show the adaptively- 
smoothed FB image and the locations of HAEs and QSOs that 
are associated with the 2QZ Clus structure. The image was 
smoothed using csmooth and was "flattened" by dividing 
the smoothed image by a smoothed exposure map. 

We constructed a Chandra source catalog by (1) search- 
ing the three standard-band images using wavdetect (in- 
cluding the appropriate 90% ECF PSF and exposure maps) 
at false-positive probability threshold of 2 x 10"^ (the equiv- 
alent detection threshold adopted by Elvis et al. 2009 in the 
C-COSMOS field); and (2) adjoining the three bandpass cata- 
logs using cross-band matching radii of 2.5 arcsec and 4.0 arc- 
sec for sources with off-axis angles of <6 arcmin and >6 ar- 
cmin, respectively. Through this procedure, we identified a 
total of 133 unique X-ray detected sources in the 2QZ Clus 
field. These 133 sources constitute our main Chandra catalog 
(see Table 1). 

We performed aperture photometry for our 133 sources in 
each of the three standard bands. For each X-ray source, we 
extracted source plus background counts 5sic and exposure- 
map values Tsic using circular apertures with radii correspond- 
ing to the 90% ECF, which were determined using the corre- 
sponding PSF-map value at the source location. We visually 
inspected these regions and found that only two source-pairs 
had very minor overlaps in the HB PSF (and no overlaps in the 
smaller SB and FB PSFs). The overlapping regions of these 
source pairs contained <1 HB counts; therefore, we did not 
perform any corrections to the photometry for these sources . 



For each source, local background counts and exposure val- 
ues were then extracted from each of the three bandpasses. 
This was achieved by creating images and exposure maps 
with all sources masked out of circular masking regions of 
radius 1.5 times the size of the 95% encircled energy fraction 
(estimated using each local PSF). Using these masked data 
products, we then extracted background counts Sbks ™d ex- 
posure values rbkc; from a larger background extraction square 
aperture that was centered on the source. The size of the 
background extraction square varied with each source and was 
chosen to contain ^ 30-100 background counts in all bands. 
For each bandpass, net source counts were computed fol- 
lowing A^= (5src-5bkg7;rc/7i,kg)/7ECF, where 7ecf is the ECF 
appropriate for the source extraction region and bandpass. We 
computed la level Poisson errors on the net counts following 
the methods described in Gehrels (1986). 

To calculate vignetting-corrected count-rates, we made 
use of count-rate maps, which were constructed by divid- 
ing the images by the exposure maps. The advantage of 
count-rate maps is that they allow for an accurate account- 
ing of the source intensity when gradients in the expo- 
sure are present in a source extraction region (e.g., near 
chip gaps, image edges, and bad pixels). We made use of 
the same source and background regions described above 
to extract on-source count-rates (/)si-c and background count- 
rates (f>bkg- Net count-rates were then computed following 
= (0src-'/'bkgAsrcMbkg)/7ECF, whcrc Asic and Abkg are the 
source and background extraction areas that contain exposure. 

For each source, we converted our vignetting-corrected 
count-rates to fluxes using conversion factors of [1.18, 0.659, 
and 1.96] xlO"" ergs cm'^ s"' (cnts s"')"' for the FB, SB, 
and HB, respectively. These factors assume a power-law SED 
with r = 1 .4 that is corrected for Galactic extinction. In Ta- 
ble 1, we provide the basic X-ray properties of the 133 main- 
catalog sources in the 2QZ Clus field. The survey reaches ulti- 
mate 5-count sensitivity Umits of «6.0 x 10"^^ ergs cm"^ s"'. 



i3.4 xlO ergs cm ^ s 



d.O X 10 ergs cm ^ s ' 



for 



the FB, SB, and HB, respectively; at z = 2.23, these limits al- 
low for the detection of a source with rest-frame 2-10 keV 
luminosity Lx ^ 10*^ ergs s"'. These limits are nicely com- 
patible with those achieved by the C-COSMOS survey, which 
reaches a factor of sal. 4-1. 7 times deeper in ultimate sensi- 
tivity (Elvis etal. 2009). 

3. RESULTS 

3.1. AGN Activity within z = 2.23 2QZ Clus Structure 

To measure the AGN activity in the 2QZ Clus structure at 
z = 2.23, we made use of the catalog of 22 HAEs from Mat- 
suda et al. (201 1) and four known z = 2.23 QSOs. After ac- 
counting for overlap in the HAE/QSO populations, there are 
a total of 23 unique z = 2.23 source candidates (see Fig. 1). 
The one QSO that does not have an HAE counterpart lies out- 
side the footprint of the narrow-band survey and is likely to 
be an HAE itself based on its selection from optical emis- 
sion lines. In the HAE survey, Matsuda et al. (201 1) selected 
candidate z = 2.23 sources that satisfied the following crite- 
ria: (1) K-UzSl > 0.251 (EWob,s > 50 A) and (2) the Ha flux 
excess significance (i.e., the ratio between H2SI excess and 
the uncertainty in the /r-H2Sl color) E > 2.5. The narrow- 
band imaging reaches a H2SI 5cr depth of Ral9.9 mag (AB). 
This selection is inherently subject to minor (at ^ 10-20% 
level) contamination from Paa emitters at z = 0.13, Pa/3 emit- 
ters at z = 0.65, and [O III]5007 emitters at z = 3.24; however. 
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TABLE 2 

X-RAY Detected HAEs and QSOs in the 2QZ Clus and C-COSMOS Fields 



Position (J2000) 




Count Rate (10"* cnts s"') 


























log/o.5-7 keV 


logLx 


log /Ha 


logi-Ha 




aj2000 


Snow 


XID 


SB 


HB 


</'2-7 kev/4'0.5-2 keV 


(ergs cm~^ s~') 


(ergs s"') 


(ergs cm~^ s"') 


(ergs s-i) 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


2QZ Clus 


10 03 23.0 


+00 07 25.0 


17 


15.1 ± 1.8 


7.3 ± 1.2 


0.48 ± 0.20 


-13.6 


44.7 






QSO 


10 03 38.3 


-1-00 18 23.8 


49 


13.0 ± 1.6 


10.9 ± 1.5 


0.84 ±0.18 


-13.5 


44.7 


-15.3 


43.3 




10 03 39.8 


-hOO 21 10.8 


58 


14.6 ± 1.8 


8.5 ± 1.4 


0.58 ± 0.21 


-13.6 


44.7 


-14.8 


43.7 


QSO 


10 03 51.6 


-hOO 15 02.1 


96 


3.0 ± 1.1 


<1.2 


>0.40 


-14.5 


43.8 


-15.2 


43.4 


QSO 


10 04 04.3 


-hOO 14 14.4 


118 


8.7 ± 1.3 


5.0 ± 1.0 


0.58 ± 0.25 


-13.8 


44.5 


-16.3 


42.3 


BzK 


10 04 09.1 


-hOO 15 45.9 


126 


2.7 ± 0.8 


<2.0 


>0.75 


-14.4 


43.9 


-16.2 


42.3 


BzK 


10 04 12.9 


-hOO 12 57.9 


129 


21.3 ± 1.9 


10.9 ± 1.4 


0.51 ±0.16 


-13.4 


44.8 


-14.7 


43.9 


QSO 


C-COSMOS 


09 58 54.6 


-1-02 14 03.6 


864 


2.3 ± 0.6 


5.2 ±0.9 


2.26 ± 0.30 


-14.0 


44.2 


-16.2 


42.4 




10 00 02.6 


+02 19 58.7 


774 


0.7 ± 0.2 


1.3 ±0.3 


1.94 ±0.41 


-14.6 


43.7 


-15.6 


43.0 


BzK 


10 00 26.6 


+01 58 23.0 


716 


0.6 ± 0.3 


1.3 ±0.4 


2.32 ± 0.56 


-14.6 


43.7 


-15.9 


42.7 


BzK 


10 00 44.2 


+02 02 06.9 


1456 


<0.5 


1.1 ±0.3 


<2.04 


-14.8 


43.5 


-16.0 


42.6 


BzK 


10 00 55.4 


+01 59 55.4 


1486 


<0.6 


0.8 ± 0.3 


<1.40 


-14.9 


43.4 


-16.1 


42.5 


BzK 


10 00 55.4 


+02 33 30.1 


1386 


2.8 ± 0.6 


<2.6 


>0.93 


-14.4 


43.9 


-15.8 


42.8 




10 00 57.5 


+02 33 45.2 


929 


2.7 ± 0.6 


<1.4 


>0.52 


-14.4 


43.8 


-15.7 


42.9 




10 01 15.2 


+01 57 13.3 


1105 


1.5 ± 0.3 


<1.9 


>1.27 


-14.7 


43.6 


-16.2 


42.4 


BzK 


10 01 29.9 


+02 17 04.1 


1261 


0.5 ± 0.2 


<1.7 


>3.21 


-15.1 


43.2 


-15.1 


43.5 




10 01 53.5 


+02 11 51.6 


451 


18.0 ± 1.4 


4.0 ± 0.7 


0.22 ±0.19 


-13.5 


44.7 


-16.0 


42.6 





NOTE. — Col.(l) and (2): Optical/near-IR right ascension (aj20oo) and declination (5j20oo): respectively. Right ascension is quoted in units of degrees, minutes, 
seconds. Declination is quoted in degrees, arcminutes, and arcseconds. Col. (3): X-ray catalog ID (XID). For sources in the 2QZ Clus field, the XID is from 
this work (see Table 1), and for the C-COSMOS sources, the XID is from Elvis etal. (2009). Col. (4) and (5): Count-rate for the SB and HB, respectively, in 
units of 10^ cnts s"' . Col. (6): Count-rate ratio between the 2-7 keV and 0.5-2 keV bands. Col. (7): Logarithm of the 0.5-7 keV source flux in ergs cm"- s"' . 
Col. (8): Logarithm of the rest-frame 2-10 keV luminosity in units of ergs s"' . Luminosities were calculated assuming each source is at z = 2.23. Col. (9) and 
(10): Logarithm of the Ha source flux (in ergs cm"^ s"') and luminosity (in ergs s"'), respectively, as provided by Matsuda etal. (2011) for the 2QZ Clus field 
and Sobral et al. (2012) for the C-COSMOS field. CoI.(l I): Notes on X-ray detected sources. 



for many cases we can make use of B-z' and z' -K colors 
(i.e., the BzK criteria) to constrain the redshifts to z ^ 1 .4-2.5 
and rule-out such contaminants. The BzK criteria includes 
{z'-K)as-{B-z')ab > -0.2 or (z'-K)p,b > 2.5 (Daddi etal. 
2004). From Matsuda etal. (201 1), 17 out of the 23 z = 2.23 
are either spectroscopically confirmed 2QZ Clus structure 
member QSOs or satisfy the BzK criteria (see Fig. 1). Interest- 
ingly, the three QSOs with B, z, K photometry do not satisfy 
the BzK selection criteria, indicating that AGN-dominated 
SEDs can be missed by these criteria. This is likely due to the 
presence of strong emission lines, and indeed the three QSOs 
lie only wO. 1-0.2 mag away from the BzK selection Une (see 
Fig. 3 of Matsuda etal. 2011). 

One HAE of the 23 z = 2.23 sources is located outside of 
the Chandra field of view. We matched the remaining 22 
sources (21 HAEs and four z = 2.23 QSOs) to our 2QZ Clus 
Chandra catalog using a 2. "5 matching radius (see annota- 
tions in Fig. 1). Given the HAE and X-ray source densi- 
ties, we estimate that this choice of matching radius will 
lead to a negligible number of random associations (i.e., a 
total of ~0.02 false-matches expected). We find a total of 
7 matches, giving an initial AGN fraction of !^2>2^y^% (la- 
intervals based on Gehrels 1986). Four of the seven X-ray 
detected sources are the optically-identified QSOs that were 
used to originally select the 2QZ Clus for study (see Matsuda 
etal. 2011 for details); therefore, three additional z = 2.23 
AGN are identified here as a result of our Chandra obser- 
vations. If we exclude the four pre-selected QSOs, the AGN 
fraction of the underlying HAE population is Six 
of the 7 X-ray detected sources are either spectroscopically 
confirmed 2QZ Clus member QSOs or BzK galaxies and are 



unlikely to be interlopers. The one HAE (XID = 49) that 
does not satisfy either criteria is inferred to be very luminous 
logLx/ergs s~' « 44.7 and may have an AGN-dominated 
SED that does not satisfy the BzK criteria (similar to the other 
known QSOs); however, we caimot rule out the possibility 
that this source is a low-redshift interloper. In Figure 2, we 
show the HB-to-SB count-rate ratio versus 2-10 keV lumi- 
nosity for X-ray-detected HAEs and z = 2.23 QSOs in the 
2QZ Clus field, and in Table 2, we Ust the properties of these 
sources. These seven X-ray sources have HB-to-SB count- 
rate ratios and X-ray luminosities that are consistent with 
powerful unobscured QSOs. 

In Figure 3, we show the distributions of Ha fluxes for the 
2QZ Clus HAEs and highlight the subset of HAEs with X-ray 
detections. We find that the AGN occupy the bright-end of the 
Ha flux distribution (/ho ~ 2 x 10"'^ ergs cm~^ s"^), indicat- 
ing (not surprisingly) that the AGN themselves have a dom- 
inant contribution to the intensity of the Ha emission lines 
in these systems. The mean Ha flux of the 2QZ Clus sam- 
ple overall appears to be higher than any of the non-AGN 
HAEs, indicating that the bulk of the Ha power from the 
2QZ Clus population is likely to be from AGN. However, 
once the pre-selected QSOs are removed, the mean Ha flux of 
the 2QZ Clus sample is within the range of the star-forming 
galaxies. A K-S test reveals that both the 2QZ Clus and C- 
COSMOS Ha flux distributions are consistent with each other 
whether or not pre-selected QSOs are included. 

3.2. Comparison with C-COSMOS 

One of our key goals is to compare the AGN activity of 
z = 2.23 HAEs in the 2QZ Clus structure with that of HAE 
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Fig. 3. — Ha flux distributions of z = 2.23 HAEs (open histograms) and 
their AGN subpopulations (filled histograms) for the 2QZ Clus (top) and 
C-COSMOS (bottom) samples. The AGN population appears to be prefer- 
entially located in HAEs with high Ha flux indicating that the HAE flux 
itself is likely influenced by the AGN. The vertical dashed line in each panel 
shows the mean Ha flux of each sample. We note that the mean Ha flux of 
2QZ Clus including pre-selected QSOs appears to be higher than any non- 
AGN HAEs, suggesting that the selection of this sample is overwhelmed 
by the AGN. When we exclude pre-selected QSOs, the mean Ha flux be- 
comes more consistent with that of the C-COSMOS field and the star-forming 
galaxy population (vertical dotted line). 

field samples. As described in §1, the 2QZ Clus structure was 
initially selected as an overdensity of QSOs, and therefore 
represents a region biased toward both active SMBH mass 
accretion and galaxy stellar growth. To put into broader con- 
text the AGN activity in the 2QZ Clus structure, we make 
use of an independent sample of z = 2.216-2.248 HAEs from 
HiZELS in the wide-area wl.6 deg^ COSMOS survey field 
(Sobral etal. 2012; see also, Geach etal. 2012). The HiZELS 
HAEs were selected using the same telescopes, cameras, and 
filters as used for the 2QZ Clus HAE selection; however, un- 
like 2QZ Clus, the COSMOS field was not selected to have 
any QSO overdensity at z = 2.23. In total there are 353 z = 2.23 
HAE candidates with /ho ^ 5 x 10~'^ ergs cm~^ s"^ (com- 
plete at the ft!90% level) that satisfy the criteria used to se- 
lect the 2QZ Clus HAEs (i.e., K-WjSl > 0.215; H2SI < 
19.9 mag). In this selection, we do not remove non-BzK 
galaxies for fair comparisons with 2QZ Clus HAEs and to 
avoid z = 2.23 AGN being removed (see discussion in § 3.1 
regarding z = 2.23 2QZ Clus QSOs that do not satisfy BzK se- 
lection). As such, our selection is somewhat more liberal than 
the final selection of galaxies presented in Sobral et al. (2012), 
which do include a BzK filtering (slightly modified), together 
with another color-color selection (UBR), high quality photo- 
metric redshifts, and information on double/triple line emit- 
ters (see Sobral etal. 2012 for details). 

Of the 353 HiZELS HAEs, 210 overlap with the footprint 
of the Chandra-COSMOS survey (C-COSMOS; Elvis etal. 
2009; Puccetti etal. 2009). C-COSMOS is a contiguous 



^0.92 deg^ Chandra survey that reaches 70-180 ks depths 
(vignetting-corrected) across s=;84% of the surveyed area. The 
combination of the HiZELS HAE and C-COSMOS surveys 
therefore constitutes a powerful data set by which direct com- 
parisons can be made with the 2QZ Clus data. We made use of 
the BzK technique to make a first-order assessment of the con- 
tamination fraction for our HAE sample. Of these 210 HAEs, 
we found that 160 had B, z, and K photometry that would al- 
low for a BzK redshift assessment. Of these 160 HAEs, 132 
(«83%) have colors satisfying the BzK color selection. We 
note that this BzK fraction is larger than that reported by So- 
bral et al. (2012) for the full HAE sample. The difference here 
is that we have limited our analysis to the relatively bright 
HAEs in the HiZELS sample, which has a lower contamina- 
tion fraction. Furthermore, the BzK fraction for our HAEs is 
consistent with that found for the 2QZ Clus HAE sample. We 
have tested the effects of excluding all non-Bz^T HAEs from 
our samples and find that this choice has an insignificant im- 
pact on the results presented throughout this paper. 

Using a matching radius of 2.5 arcsec, we obtain a total of 
10 Chandra sources matched to the 210 HiZELS HAEs that 
were within the C-COSMOS footprint, giving an initial AGN 



fraction of 4.8 



+2.0c 
■1.4 ' 



a factor of 6.7;^ 2 times lower than the 



initial 2QZ Clus AGN fraction of 32!}^%. Given the QSO pre- 
selection in the 2QZ Clus field, this difference in AGN frac- 
tion may not be so surprising. When we exclude the four pre- 
selected QSOs from our computation, we find the 2QZ Clus 
AGN fraction is « 17!^''% (i.e., 3 AGN out of 18 HAEs); still a 



factor of 3. 5^2 1 times higher than the HiZELS AGN fraction. 
If we assume that the C-COSMOS HAE AGN fraction is rep- 
resentative of the HAE population in general, we estimate that 
the binomial probabiHty of detecting 3 or more AGN out of 
18 HAEs is «6%. 

In Figure 2, we show the count-rate ratio vs. 2-10 keV lu- 
minosity for the HiZELS z = 2.23 HAEs that are X-ray de- 
tected in C-COSMOS. The HiZELS HAEs host an almost 
equal blend of moderately obscured (A^h ~ 10^^ cm^) and un- 
obscured AGN as measured by their band ratios. This ap- 
pears to differ somewhat from the "softer" AGN found in 
2QZ Clus (see Fig. 2). If we assume that the ~50% ratio 
of obscured to unobscured AGN in the C-COSMOS HAEs is 
typical of z = 2.23 HAEs, then the probability of finding aU 
three uniquely X-ray selected AGN (i.e., excluding the four 
pre-selected QSOs) in 2QZ Clus to be unobscured is w25%. 
We therefore do not conclude that the differences between the 
X-ray spectral slopes in 2QZ Clus HAEs are significantly dif- 
ferent from those of C-COSMOS HAEs. In Figure 3, we 
display the Ha flux distribution of the C-COSMOS HAEs 
and highlight the subset that are X-ray AGN. In contrast to 
2QZ Clus HAEs, we find that the mean C-COSMOS HAE 
flux is within the distribution of non-AGN HAEs, indicating 
that the total Ha power output is likely dominated by the star- 
forming galaxy population. 



4. DISCUSSION 

The above analyses indicate that the 2QZ Clus HAEs con- 
tain an enhanced AGN fraction compared to HAEs found 
in the C-COSMOS field. Here we discuss the possible cir- 
cumstances surrounding the enhanced SMBH growth in the 
2QZ Clus structure. We further put into context the mu- 
tual SMBH and galaxy growth of 2QZ Clus and C-COSMOS 
HAEs. 



Concurrent SMBH and Galaxy Growth in z = 2.23 HAEs 
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TABLE 3 
X-RAY Stacking of HAE Samples 



logLx Mbh S250^im SFR 

Sample \og(p/{p)) N.,,,, Nmk ,/agn (^-t kov A\t 5-2 kev (ergss"') (ir^Mo yr') (tnJy) (Mo yr') 

(1) (2) (3) (4) (5) (6) (7) (S) (9) (10) 

C-COSMOS Low Density -0.38 ± 0.38 61 2 0.03;j| "^ 0.62 ± 0.47 42.52 ±0.13 11 ±4 3.9 ± 1.6 68 ± 29 

Medium Density 0.38 ± 0.38 122 5 0.04-^%l 0.53 ± 0.66 42.79 ± 0.23 21 ± 15 3.4 ± 0.6 60 ± 10 

High Density 1.12 ±0.38 27 3 O.ll+Jj,^^ 0.45 ± 1.03 42.50 ±0.17 11 ±5 <2.3 <42 

2QZ Clus including QSOs 0.31 ±0.10 21 6 0.29!^-J] 0.48 ±0.10 44.01 ±0.14 362 ± 142 6.0 ± 2.6 111 ±73 

excluding QSOs 0.31 ±0.10 18 3 0-17!^:nq 0.59 ±0.18 43.76 ± 0.19 204 ±110 3.6 ± 1.3 63 ±23 



NOTE. — Col.(l): Description of z = 2.23 HAE sample being stacked. Col.(2): HAE source overdensity, p/ (p), computed as the ratio of the local HAE source 
density p to mean HAE density across the entire HiZELS COSMOS survey area (p) (Sobral et al. 2012). For the shallower 2QZ Clus HAE sample, we computed 
p/ (p) using bright HAEs in both the 2QZ Clus and HiZELS COSMOS fields (see §4 for details). Col. (3): Number of galaxies in the relevant stacking sample. 
This number reflects exclusion of HAEs that were outside the C-COSMOS footprint, as well as sources that were in the near vicinity (within ail5 arcsec) of 
unrelated X-ray detected sources. Col.(4): Number of sources that were detected in the X-ray band. Due to the relatively high survey luminosity limits for 
z = 2.23, these sources are expected to be AGN. Col.(5): AGN fraction of sample. Col. (6): Ratio of stacked 2-7 keV to 0.5-2 keV count-rates <f). Col. (7): 
Logarithm of the mean 2-10 keV luminosity calculated using the stacked 0.5-2 keV emission. Col. (8): Estimate of the mean BH accretion rate based on the 
X-ray luminosity (see §4 for details and assumptions). Col. (9): Mean 250/im flux in mjy based on stacking analyses. Col. (10): Mean SFR based on 250/im flux 
(see §4.2). 
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Fig. 4. — (a) AGN fraction versus local source overdensity for the C-COSMOS (filled circles) and 2QZ Clus (filled squares) z = 2.23 HAEs. For the 2QZ Clus 
field, we have highlighted results for both the inclusion (black squares) and exclusion (gray squares) of the QSOs that were identified via the 2dF QSO Redshift 
Survey and used to select the 2QZ Clus field. We find suggestive evidence for an increase in AGN fraction with local density for the C-COSMOS HAEs and find 
that the 2QZ Clus AGN fraction is significantly elevated over that of the C-COSMOS HAEs. (b) Mean SMBH mass accretion rate per SFR vs. local overdensity 
for each of the stacked samples presented in § 4.2 and Table 3. For comparison, the local MBu/^gal relation and its dispersion from Haring & Rix (2004) has 

been shown as a dotted line with shaded region. The Mbh/SFR values for 2QZ Clus HAEs appear to be consistent with the local MsH/A^gai relation, while the 
C-COSMOS HAEs fall below this relation. We are unable to determine whether Mbh/SFR varies with local environment for the C-COSMOS HAEs. 



4. 1. Environment as a Driver of Galaxy and SMBH Growth 

To assess the role that environment plays in driving SMBH 
accretion in the HAE population, we made use of the wide- 
area HiZELS HAEs in the C-COSMOS survey to estimate 
AGN fraction as a function of source density. For the broader 
HiZELS HAE sample (i.e., the « 1.6 deg^ sample), we mea- 
sured local HAE overdensities for each source following 
p/{p) ~ A / (jxrl) / {p) , where is the separation between a 
given source and its fourth nearest neighbor and {p) is the 
mean HAE source density over the entire HiZELS COS- 
MOS field. Since the HiZELS COSMOS field has wider 
aerial coverage than the subpopulation of HAEs that lie within 
the C-COSMOS footprint, we expect that p/{p) should pro- 
vide a good estimate of the local source environment for all 
C-COSMOS sources without suffering from uncertainties due 
to survey edge effects. Focusing on only the HAEs that were 
within the C-COSMOS footprint, we divided the sample into 



three bins of source overdensity, log p/{p) = [-0.35, 0.35, 
1 . 12] ±0.38. For each bin, we computed the fraction of HAEs 
hosting luminous AGN detected in the X-ray band. Figure 4a 
shows the AGN fraction as a function of overdensity p/{p) for 
the C-COSMOS sources, and in Table 4 we tabulate the AGN 
fractions of each subset. We find suggestive evidence for an 
increase in AGN fraction with local HAE source density, such 
that the /AGN(logp/(p) « 1.12) = 3.4!^:° x /AGN(log/o/(p) « 
-0.4). A similar trend was noted for Lyman- a emitters in the 
z = 3.1 SSA22 protocluster (Lehmer etal. 2009a). 

Given that the 2QZ Clus HAE survey is somewhat shal- 
lower and less complete than the HiZELS COSMOS HAE 
survey, we cannot directly compute the equivalent HAE 
source density all the way down to the Ha flux Umit. There- 
fore, to estimate an average p/ (p) for all sources in 2QZ Clus, 
we computed the survey-wide 2QZ Clus HAE source den- 
sity of a highly complete ( ^ 90%) subsample using only 
relatively bright HAEs (fua > 6.3 x lO"''' ergs cm"^ s"'). 



8 



LEHMERETAL. 



and then compared this with the source density of relatively 
bright COSMOS HiZELS HAEs. The relative bright-source 
densities indicate that the 2QZ Clus HAEs have on average 

(P)2QZ cius/(/'>h"Iels = 2-0 ± 0.5. This value is somewhat 
higher than, although consistent with, the «1.5±0.4 value 
found by Matsuda et al. (201 1) when comparing the 2QZ Clus 
HAEs with a smaller HAE control sample. We note that the 
relatively small area in 2QZ Clus limits us from effectively 
measuring the local densities of each galaxy in the way we 
did for the HiZELS HAEs. A wider area HAE survey of the 
2QZ Clus structure would mitigate this limitation. In Fig- 
ure 4a, we highUght the AGN fraction of the 2QZ Clus HAEs, 
both including and excluding HAEs that were pre-selected to 
be QSOs. After excluding the pre-selected QSOs, we find 
that the 2QZ Clus AGN fraction is a factor of w4.1+2 ° times 
higher than the AGN fraction of C-COSMOS HAEs in simi- 
larly overdense environments. 

From the above analyses, it appears as though the enhanced 
AGN fraction in the 2QZ Clus cannot be explained by an en- 
hanced source density. However, given that the HAE source 
selection is limited to the detection of sources with high-SFR 
(SFR >, 14 Mq yr"') or AGN activity, it may be possible that 
there are large numbers of 2QZ Clus galaxies that are lower 
SFR and/or high Ha extinction leading to an inaccurate char- 
acterization of the local source environment. For this to ex- 
plain the relatively large AGN fraction in terms of environ- 
ment would require that the underlying galaxy population in 
2QZ Clus be very different from that of C-COSMOS. Future 
observations sensitive to these galaxy populations would be 
needed to address this. 

If we assume that the enhanced AGN fraction in the 
2QZ Clus structure is due solely to the biased QSO-based 
selection and that our estimates of the local source overden- 
sities are representative of the true environments in which 
the galaxies are found, we can assess how rarely we should 
find such elevated AGN fractions. From our analysis of the 
C-COSMOS HAEs, we know that the probability of an HAE 
hosting an AGN in an environment like that of 2QZ Clus is 
s=;4%. Assuming this incidence fraction is representative of 
the broader HAE population, we used the binomial probabil- 
ity distribution to compute the rarity of finding AGN frac- 
tion enhancements like those of 2QZ Clus. We find that the 
probability of detecting 3 or more AGN in 18 HAEs (i.e., af- 
ter removing pre-selected QSOs) is «4.1%. This suggests 
that there is likely to be something inherently "special" about 
the 2QZ Clus region; perhaps some unaccounted for physical 
mechanism is responsible for the enhanced AGN activity that 
is unique to 2QZ Clus. It is also possible that the elevated 
AGN fraction is a product of the presence of more massive 
host galaxies and SMBHs. For example, Xue etal. (2010) 
found that the AGN fraction increases rapidly with stellar 
mass. From their AGN fraction versus stellar mass relations 
for luminous AGN (see Fig. 14B of Xue etal. 2010), we in- 
fer that the 2QZ Clus galaxies would have to be on-average 
a factor of w5-8 times more massive than the C-COSMOS 
galaxy population to explain the enhanced AGN fraction that 
is observed. The 2QZ Clus and C-COSMOS HAEs were se- 
lected at the same SFR limits and have similar mean SFRs 
(see below). Since SFR is in general correlated with stellar 
mass, with a dispersion of a factor of «2 (e.g., Elbaz etal. 
2007; Salim etal. 2007), it is unlikely that the mean stellar 
masses of the HAEs in 2QZ Clus would be «5-8 times higher 
on-average than those in C-COSMOS. Better characterization 



of the 2QZ Clus HAE source properties, large-scale environ- 
ment, and underlying non-HAE galaxy population will help 
us understand this. This could be achieved via HST WFC3 
infrared observations of the HAE population and/or spectro- 
scopic follow-up of BX/BM and BzK galaxy candidates. 

4.2. Relative SMBH and Galaxy Growth Rates 

To estimate the mean SMBH growth rates of 2QZ Clus 
and C-COSMOS HAEs, we performed X-ray stacking fol- 
lowing the techniques described in Lehmer et al. (2008), with 
errors measured as 1 cr intervals of bootstrap resampling (see, 
e.g., Basu-Zych etal. 2012 for details). In 2QZ Clus, we 
stacked the total HAE samples both including and exclud- 
ing the pre-selected QSOs. For C-COSMOS, we stacked 
each of the three subsamples that were divided based on 
source density (see §4.1). In Table 3, we tabulate our stack- 
ing results for the five subsamples. We find that the stacked 
2-10 keV luminosities range from logLx/(ergs s"') « 42.5- 
44, indicating that AGN clearly dominate the mean stacked 
X-ray emission over normal galaxy emission by more than 
an order of magnitude for galaxies of similar SFRs and red- 
shifts (e.g., Lehmer etal. 2010; Basu-Zych etal. 2012). In 
all cases the 2-7 keV/0.5-2 keV count-rate ratios indicate the 
mean stacked spectrum is consistent with unobscured AGN 
(see Fig. 2 and Col. 6 in Table 3). We converted our mean 
2-10 keV luminosities into bolometric luminosities L^^^ us- 
ing a correction factor of 22.4, which corresponds to the me- 
dian bolometric correction for local AGN with Lx ~ 10"*^- 
lO'*^ ergs s~' (Vasudevan & Fabian 2007). We then esti- 
mated the mean SMBH accretion rates following Mbh ~ 
(1 - e)L^^^ /(ec-), where c is the speed of light and e is the 
efficiency by which accreting mass is converted into radia- 
tive energy; here we assume e ~ 0.1 (see, e.g., Marconi etal. 
2004 for motivation). These assumptions yield mean mass 
accretion rates spanning Mbh ~ 0.01-0.02 Mq yr~' for the 
C-COSMOS HAE samples and «0.2-0.4 M© yr"' for the 
2QZ Clus HAE samples (see Col. 8 in Table 3). 

The above X-ray stacking shows that the mean SMBH 
accretion rates of 2QZ Clus HAEs are significantly larger 
than those of C-COSMOS HAEs; however, we want to test 
whether the corresponding SFRs of the HAEs are on average 
consistent with that expected from the local MBn/Mgai rela- 
tion. Although the Ha emission line provides an immediate 
proxy for the star-formation rate (SFR), Ha is often subject 
to extinction in starburst galaxies like those studied here (e.g., 
Calzetti etal. 1997). Furthermore, in the case of QSOs, the 
Ha emission line fluxes from the AGN will overwhelm any 
signal related to star-formation activity, making it difficult to 
infer directly the SFRs of these sources (see, e.g.. Fig. 3). 

To remedy these issues, we have made use of Herschel 
observations at 250/um using the Spectral and Photometric 
Imaging Receiver (SPIRE). Herschel SPIRE data are avail- 
able in the C-COSMOS field via the HerMES campaign 
(Oliver etal. 2012) and in the 2QZ Clus field through a GO 
program (PI: Matsuda; Matsuda etal. 2013, in-preparation). 
These data allow us to probe emission from the rest-frame 
80 fxm, which is well beyond the peak of the expected dust- 
emission associated with QSOs (SS6-15 ^m; e.g., Efstathiou 
& Rowan-Robinson 1995; Netzer et al. 2007; Hatziminaoglou 
et al. 2010; Mullaney et al. 201 1). We therefore anticipate that 
SPIRE probes the cool dust emission associated with star- 
formation activity, and the total SFR can be measured us- 
ing these data. We performed SPIRE stacking of the three 
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Fig. 5 . — Mean SMBH mass accretion rate per SFR (Mbh/SFR) versus 2-10 keV luminosity for the stacked values of 2QZ Clus and C-COSMOS HAEs (see 
Table 3). In the left panel, we compai'e these values with z~2 AGN-selected samples from Hanison etal. (2012; H12) and z ~ 1-3 AGN from Mullaney etal. 
(2012a). In the right panel, we compare 2QZ Clus and C-COSMOS A/bh/SFR values with those of z~ 1-2 star-forming galaxies from Mullaney etal. (2012b) 
and the range of vigorously star-forming submilHmeter galaxies from Alexander etal. (2008). For reference, we have shown in both panels the local MgH/^eal 
relation and its dispersion as repoiled by Haring & Rix (2004). The 2QZ Clus HAEs have mean Mbh/SFR values consistent with the trends observed for AGN, 
while the C-COSMOS HAEs appear to be broadly consistent with both the z ~ 1-2 AGN and star-forming galaxy populations. 



C-COSMOS and two 2QZ Clus samples described in Table 3 
following the techniques detailed in Harrison et al. (2012). We 
obtained significant signal in the 250/im stacks for all sam- 
ples except for the highest overdensity C-COSMOS subsam- 
ple (see Table 3). Errors on mean SPIRE flux represent the 
largest Icr interval based on bootstrap resampling (see Har- 
rison etal. 2012 for details). We converted 250/im fluxes to 
LiR (8-1000 ^m) using the SED library of Chary & Elbaz 
(2001), selecting a SED (redshifted to z = 2.23) on the basis 
of the total monochromatic luminosity probed by the 250/im 
emission. We converted the infrared luminosities to SFRs fol- 
lowing SFR/(Mq yr"') w 9.8 x 10"'' x Lir/Lq. This conver- 
sion assumes negligible contributions from UV emission and 
is appropriate for a Kroupa (2001) IMF (see eqn. 1 of Bell 
etal. 2005). The mean SFRs are reported in Column 10 of 
Table 3. 

From the above X-ray and infrared stacking analyses, 
we obtain mean Mbh/SFR ratios for our samples. In 
Fig Ab, we display Mbh/SFR versus overdensity {p/ (p)) for 
these samples. We find that the 2QZ Clus HAEs have 
Mbh/SFR « 3.2 x lO"-', similar to the mean MBn/Mgai ratio at 
z = (e.g., Haring & Rix 2004; Bennert et al. 201 1). For the 
C-COSMOS HAE subsamples, we find much lower values 
of Mbh/SFR w (0.2-0.4) xlO'l Such a deficit in SMBH- 
to-galaxy growth rate ratio may be due to the SFR-biased 
selection of C-COSMOS HAEs, which has constrained our 
analysis to include only powerful star-forming galaxies with 
extinction uncorrected Ha-based SFRhq > 14 Mq yr"'. 

Recent Herschel SPIRE analyses of X-ray selected AGN 
in the CDF-N, CDF-S, and COSMOS fields have shown that 
on average the mean SFR increases with mean AGN lumi- 
nosity; however, there is some debate as to whether the most 
powerful QSOs with Lx ^ lO"*^ ergs s"' have declining SFRs 
(e.g., Harrison etal. 2012; Page etal. 2012). In the left panel 
of Figure 5, we compare the 2QZ Clus and C-COSMOS 
HAE Mbh/SFR values with those of X-ray selected AGN 
presented in Harrison etal. (2012) and z ~ 1-3 AGN from 
Mullaney etal. (2012a). X-ray luminosities have been con- 



verted to Mbh following the techniques discussed above and 
SFRs have been converted to be consistent with our choice 
of Kroupa (2001) IMF. It is appai-ent that the mean Mbh/SFR 
for X-ray selected AGN appears to be consistent with the lo- 
cal MBn/Mgai relation (Haiing & Rix 2004) for AGN with 
Lx ^ lO"'^ ergs s~'. At higher Lx, MBn/Mgai increases dra- 
matically. The 2QZ Clus HAE samples have Mbh/SFR vs. Lx 
mean values in good agreement with the trends observed for 
the Harrison etal. (2012) X-ray selected AGN (see Fig. 5). 
As discussed above, the C-COSMOS HAEs appear to have 
somewhat lower Mbh/SFR values below that expected from 
the local MBn/Mgai relation; however, these values appear to 
be consistent with those of the low-luminosity X-ray selected 
AGN from Mullaney et al. (2012a). Given that the mean SFRs 
and masses of all the AGN samples shown in the left-panel 
of Figure 5 are relatively high (SFR w 10-100 M© yr"' and 
M,t ^ 5 x 10'°), this trend suggests that SMBHs and galaxies 
do not grow simultaneously. 

In the right-panel of Figure 5, we compare the 2QZ Clus 
and C-COSMOS HAE Mbh/SFR values with those of z « 1-2 
star-forming galaxies (Mullaney etal. 2012b). We find that 
the C-COSMOS HAEs have Mbh/SFR in good agreement 
with field z ~ 1-2 star-forming galaxies, which appear to also 
agree with z ~ 2 submilHmeter galaxies studied by Alexander 
etal. (2008). These values appear to lie a factor of sa3 below 
the local MBn/Mgai. We note, that the various assumptions 
(e.g., bolometric correction and accretion efficiency) and un- 
known uncertainties used in our calculations may contribute 
to systematic offsets in Mbh/SFR. However, if our findings 
are correct, typical star-forming galaxies must undergo vig- 
orous SMBH growth in short duty cycles, where Mbh/SFR 
is above the local MBn/Mgai ratio. If the QSO phase at 
Lx ^ 10"'^ ergs s~' is responsible for such episodic growth, 
we can use our observations to estimate the QSO duty cycle 
required to bring star-forming galaxies onto the MBn/Mgai re- 
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lation: 

Mbh/M^,i « (MBH/SFR)Qso/Qso + (MBH/SFR)gai(l -/qso), 

(1) 

where /qso is the fraction of time that a galaxy is in a QSO 
phase. Taking Msh/M^,i = IQ-^ ^-lQ-^ o, (Mbh/SFR)qso « 
IQ-i ^ and (MBH/SFR)gai « lO'^-^, we find /qso « 2-8%. For 
C-COSMOS HAEs, we estimate that the QSO fraction for 
Lx ^ 10^ ergs s'^ AGN is «1.0!i|%. This value and our 
duty-cycle estimate of /qso are consistent within errors, indi- 
cating that the short-term SMBH growth in a QSO phase may 
be sufficient to bring the typical HAE up to the Msu/Mgai 
relation. This is consistent with the observation that the ma- 
jority of SMBH growth density (using X-ray emissivity as a 
proxy) in the Universe at z 2 can be attributed to luminous 
AGN (Lx ^ 10'*'* ergs s"'; e.g., Hasinger etal. 2005); a sig- 
nificant fraction of this population is expected to be fueled by 
mergers (~30^0%; e.g., Treister etal. 2012). We note that 
none of the luminous AGN among the C-COSMOS HAEs 
fell within the observational footprint of the HST WFC3 IR- 
imaged region of the the Cosmic Assembly Near-infrared 
Deep Extragalactic Legacy Survey (CANDELS; Grogin etal. 
2011; Koekemoer etal. 2011), so it was not possible to ex- 
amine directly the rest-frame optical morphologies of the 
HAE AGN subpopulation. We used the HST ACS optical 
(F814W) morphology catalogs from Tasca etal. (2009) from 
the broader COSMOS field (Scoville etal. 2007) to com- 
pare the concentrations and asymmetries of the rest-frame UV 
emission from the HAEs and their AGN subpopulation. These 
parameters were available for 270 C-COSMOS HAEs and all 
10 AGN. Compared with the C-COSMOS HAE population, 
the AGN subpopulation morphology distribution is skewed 
towards higher optical-light concentrations (K-S test reveals 
the populations differ at the 98.1% confidence level). The 
AGN subpopulation asymmetry distribution is consistent with 
the broader HAE population, with no evidence for more lumi- 
nous AGN having larger asymmetries, as might be expected 
from merging systems. However, the observed HST images 
reveal that the rest-frame UV light from the most luminous 
AGN is dominated by the QSO emission, skewing the light to 
higher concentrations and making merger signatures elusive. 
Observations of these sources with HST WFC3 IR (rest-frame 
optical) would improve our ability to address whether merg- 
ers play a significant role in triggering driving SMBH growth 
in the HAE population. 

5. SUMMARY 

We have conducted a «100 ks Chandra observation over 
the 2QZ Clus structure and present a Chandra point-source 
catalog. Within the Chandra footprint, the 2QZ Clus struc- 
ture contains 21 HAEs and four QSOs at z = 2.23 (total of 
22 unique z = 2.23 sources) and is on-average a factor of «2 
times overdense compared to field HAEs. The 2QZ Clus 
was initially selected by (1) identifying an overdensity of 
Z K, 2.23 QSOs and (2) performing Ha narrow-band photom- 
etry to preferentially detect additional star-forming galaxies 
and AGN at z « 2.23 (via HAE selection; see Matsuda etal. 
2011). Given its selection, the 2QZ Clus HAEs contain a 
rich mixture of AGN and star-forming galaxies. To put into 
broader context the 2QZ Clus HAEs, we compare their prop- 
erties with those of a larger sample of 210 z = 2.23 HAEs se- 
lected in the C-COSMOS field. Our findings are summarized 
below. 



• We find seven of the 22 z = 2.23 sources are detected in 
the X-ray band, including all four QSOs. These sources 
have 2-10 keV luminosities of «(8-60) x lO'*^ ergs s~' 
and HB/SB count-rate ratios indicative of unobscured 
QSOs. Comparison with HAEs in C-COSMOS reveal 
that the HAE AGN fraction in the 2QZ Clus field is 
enhanced by a factor of Ki'i.5^2 \ compared with the 
broader field (after excluding pre-selected 2QZ Clus 
QSOs). Therefore, the 2QZ Clus structure is a region in 
which SMBH growth is in a particularly active phase. 

• On average. X-ray selected AGN and QSOs occupy the 
high-flux end of the Ha flux distribution, suggesting 
that the Ha line for these bright HAEs is on average 
dominated by the QSO component. For 2QZ Clus the 
mean Ha flux is higher than any non-AGN HAE, in- 
dicating that the total 2QZ Clus Ha power output for 
HAEs is dominated by AGN/QSO activity. However, 
the C-COSMOS HAEs have Ha power dominated by 
star-forming galaxies. 

• We find suggestive evidence that the C-COSMOS 
z = 2.23 HAE AGN fraction increases with increas- 
ing local HAE overdensity. The 2QZ Clus HAEs re- 
side in a factor of «2 HAE overdense environment. 
After excluding known optically-selected QSOs, the 
2QZ Clus AGN fraction is a factor of «4.1+f;g times 
higher than C-COSMOS HAEs in similarly overdense 
enviroimients. Either the local density as measured by 
HAEs alone does not provide an accurate characteriza- 
tion of a truly more overdense 2QZ Clus environment, 
or the 2QZ Clus HAE population is a statistical out- 
lier (at the «95.9% level) in terms of its elevated AGN 
fraction. 

• We make use of Chandra and Herschel 250/im stack- 
ing of 2QZ Clus and C-COSMOS HAE samples to 
measure mean SMBH mass accretion rates and SFRs. 
These stacking analyses reveal that the mean X-ray 
emission from all HAE samples is dominated by AGN, 
with 2QZ Clus HAEs having QSO-hke luminosities 
(i-a-iokev « [6-10] xl043 ergs s'^) and C-COSMOS 
HAEs having much lower Seyfert-like luminosities 
(^2-10 keV ~ [3-6] xlO'*^ ergs s"'). The mean SFRs 
range from «60-110 M© yr^^ for 2QZ Clus HAEs 
depending on whether optically-selected QSOs are re- 
moved or not. The SFRs for C-COSMOS HAEs span 
^40 Mq yr~^ to «70 Mq yr~'. The inferred mean 
SMBH mass accretion rate to SFR ratio (Mbh/SFR) 
for these samples indicate that the 2QZ Clus SMBHs 
and galaxies appear to be growing at rates comparable 
to those expected by the local Msn/Mg^i ratio, while 
the mean Mbh/SFR for C-COSMOS HAEs he a factor 
of w3 times lower than the local M^n/M^^x ratio. We 
find that the Mbh/SFR vs. X-ray luminosity for both the 
2QZ Clus and C-COSMOS HAE samples appear to fol- 
low trends found for z « 1-2 star-forming galaxies and 
z i=a 2 X-ray selected AGN. 

• We estimate that an episodic QSO phase with a duty 
cycle of «2-8% would allow C-COSMOS HAEs to 
emerge onto the Msn/Mg^i relation despite their lower 
population-averaged Mbh/SFR. This estimate is consis- 
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tent with the observed C-COSMOS HAE AGN fraction 
(wO.4-2.3%) for Lx ^ 10'*'* ergs s"' sources. 



sive passive galaxies and additional galaxies with low SFRs 
or heavily extinguished Ha emission). 



Future observations of the 2QZ Clus are needed to further 
evaluate (1) what physical conditions have led to such an en- 
hanced AGN fraction, (2) the large-scale environment of the 
2QZ Clus structure, and (3) whether there is a significant un- 
derlying overdensity of galaxies that are not HAEs (e.g., mas- 
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